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PREFACE 

WHILE lecturing in 1913 on the calculation of self and mutual 
induction, the writer noticed that the integral in the case of co- 
axial circles was of such a form as to make the application of the 
arithmetico-geometrical means obvious. An extremely elegant formula 
(Proc. Roy. Soc. A, Vol. c. 1921, p. 63. Appendix, this book, No. 35) 
when put to numerical test justified the method as one well adapted to 
computation. The calculation of more complex cases of self and mutual 
induction led to several interesting formulae for the complete third 
elliptic integral. With a view to making a complete study of the group 
of recurrence formulae associated with the A.O.M. scales, further work 
was postponed until a vacation in 1921, when many of the results in- 
cluded in this book were obtained, and were thought by the 'writer to 
be new. On later consulting the Collected Works of Legendre, Gauss 
and Jacobi, the writer found that he had rediscovered many known 
formulae, but that quite a number, likely to be of service in computation, 
appeared to be hitherto unknown. In the summer of 1922 the writer 
decided, as a vacation task, to include in a single monograph the entire 
theory of elliptic functions associated with the A.G.M. scales, thus com- 
pleting, and to some extent adding to, the work of Gauss left unpublished 
after his death, and placing before mathematicians in accessible form 
a mode of approach to Elliptic Function Theory directly related to the 
art of machine-computation. It is interesting to note in this connection 
that this subject was to have formed the content of the third volume 
of Halphen's Fonctions Ettiptiqws, unfortunately left incomplete at the 
time of his death. ' 

While the present monograph contains many new formulae and 
methods of computation, no claim is made as to novelty in fundamental 
analytical treatment. References to more important books and memoirs 
consulted by the author are given in foot-notes. A complete bibliography 
is not given since the reader will find exhaustive references in the Royal 
Society Index, 1800-1900, Vol. I, Mathematics, and after 1900 in tha 
mathematical volumes of the International Catalogue ~qf Sciential 
Literature. 
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ON THE DIRECT NUMERICAL CALCULATION 
OF ELLIPTIC -FUNCTIONS AND INTEGRALS 

SECTION I 
INTRODUCTION 

ALTHOUGH the numerical evaluation of elliptic functions and integrals 
is dealt with in many of the standard treatises, it cannot be said that 
the use of tables is entirely satisfactory in dealing with computations 
required in several branches of physics and astronomy. Many of the 
integrals appearing in these problems can only be expressed as some- 
what complicated expressions involving the complete elliptic integrals 
of the first and second kinds, and the corresponding functions of com- 
plementary modulus. These formulae give rise in many cases to 
differences of nearly equal quantities : the numerical calculation of such 
expressions to a sufficient number of significant figures is thus often a 
matter of some difficulty, requiring tedious and sometimes uncertain 
interpolations. 

While the numerical calculation of the Theta-functions is extremely 
rapid by the use of the highly convergent series involving the nome </, 
these are not usually the functions which make their appearance directly 
in physical or astronomical problems. The necessary transformations 
required to express the results in terms of the ^-series require in many 
cases somewhat complicated analysis, while numerical computation by 
this method necessitates the exact calculation of a largo number of 
auxiliary quantities*. 

The forthcoming tables in course of construction under the auspices 
of the British Association t are so arranged as to allow a large number 
of functions to be evaluated, and will very greatly aid numerical work 
in many practical problems involving the use of elliptic integrals and 
functions. While sufficient for most purposes, the use of these tables 
will still involve considerable interpolation when it is required to take 
out values corresponding to arguments and moduli expressed by six or 

* On this point see Rosa, E. H., and Grover, F. N., Bulletin 169, Bureau of Stan- 
dards, Washington (1913), pp. 67 and 73; Nagaoka, H., Tok. Coll. tfc. J. 27 (1909), 
No. 6, and Jour. Tok. Math. Phys. Soc. 4 (1908) , p. 284 ; 6 (1911), p. 10; Nagaoka and 
Sakurai, Scientific Papers of the Institute of Physical and Chemical Research, Vol. n. 
Dec. 1922, pp. 1-67; Olshausen, G. R., Phys. Rev. 51 (Dec. 1910), pp. 617-636. 

t See report (by Sir A. G. Greenhill), Brit. Ass. Report, 1912 (Dundee), pp. 
39-55. 
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seven significant figures, as is now necessary in several modern applica- 
tions of elliptic functions. 

For several reasons it is advantageous to have available a direct 
method of computation, independent of auxiliary tables, based on the 
use of the modern high-speed calculating machine. The writer has 
therefore thought it desirable to set out in detail the method of 
calculation based on Landen's quadric transformation and the use of 
the scale of arithmetico-geometrical means, recapitulating briefly the 
original developments of Lagrange, Legendre, Gauss, Jacob! and others, 
extending them in several directions and keeping in mind the ultimate 
use of the calculating machine for numerical work. 

SECTION II 

HISTORICAL NOTE ON LANDEN'S TRANSFORMATION AND 
THE VARIOUS SCALES OP MODULI AND AMPLITUDES 

The importance of realizing rapid and accurate methods of calculating 
the elliptic integrals now denoted by 

it = F(<t>, k) = Jf d*/A (*, t) and E(+, k] = A <>, 4) d<f> . . .(1), 

where A (<, *) = J( I - A sin 2 </>), 

was first remarked by Euler* (1766), although it was not until several 
years later that Laudenf (1775) discovered in geometrical form the 
transformation which forms the basis of the existing methods of the 
numerical calculation of the elliptic integrals. 

A method of successive transformations for the ultimate reduction of 
the algebraic forms of these integrals to elementary integrals was pub- 
lished by Lagrange I in 1784-85. This memoir contains an exposition of 
the scales of arithmetico-geometrical means, together with two types of 
algebraic transformations corresponding to the increasing and diminishing 
amplitudes of Landen's trigonometrical forms: it is shown that the 
given integral must lie between two limits which may be made to 
approach each other as closely as may be desired : the limits thus ob- 

* Euler, Nor i Comm. Acad. Sc. Petrop. vol. x. 17(56. 

f Landen, Phil. Trans. Roy. Soc. vol. LXY. p. 283, 1775; Mathematical Memoirs, 
London, 1780. 

Lagrange, "Sur une nouvelle m&hode de calcul integral pour lea diffdrentielles 
ftltectces d'un radical carre" sous lequel la variable ne passe pas le quatrieme degrc* t " 
jl/?m. de I'Acad. roy. de Sc. de Turin 1 1. n. 1784-5; CEuvie* (Gauthier-Villars, Paris, 
1868), t. ii. pp. 253-312. 



NOTE ON LANDEN S TRANSFORMATION 3 

tained for the integrals of the first and second kinds are practically 
identical with the limits and series obtained about the same time by 
Legendre. 

Landen's transformation was applied by Legendre* to the numerical 
calculation of the elliptic integrals, and furnished the method of com- 
putation by means of which the latter's well-known tables were con- 
structed. The scales of moduli and amplitudes employed by Legendre 
are briefly described below. 

If we form successive moduli and amplitudes according to the 
recurrence formulae 

* + i -(I- *.')/(! +*) and tan (*.,- *,) = *,' tan*. ...(2), 
Landen's transformation in trigonometrical form leads to the result 

J?(**.,*.+0 = (i+*.')*'(*.,*0 .............. (3). 

As n increases, the successive moduli k n converge to zero, and the 
complementary moduli kn to unity, while the amplitudes increase so 
that < n /2 n tends rapidly to a definite limit : ultimately we may write 
A 7 (< n , kn) - <t> n and thus obtain the result 

F(+*>**y-*&KI*)Wr, where AT=Jir(l +^(1 +*,) - (1+*.) 

......... (4). 

The result of calculating successive moduli and amplitudes in reverse 
order is equivalent to calculating in forward order these quantities 
according to the recurrence formulae 

*,,+! = 2 V*/(l + *0 and sin (2i/vn - ^) = k n sin f n ...(5). 
This transformation then gives 



As n increases the successive moduli k n+l converge to uuity and the 
amplitudes to a limit i/v. we thus obtain, when n is sufficiently large, 

/'(*,*.)= I "sec <M<A B = log tan (Ja- +</,,,) ...... (7), 

/O 

giving finally 

/^(ife, *b) = */(*!*.* -. k n )Uk Q . log tan QTT + *M .-(8). 
Legendre also developed for purposes of numerical calculation extremely 
convergent series for E(<t>, k) which he employed in the tabulation of 
this function. 

* Legendre, "M^moire sur les integrations par arcs d'ellipse" and "Second 
mchnoire, etc.," Mm. de VAcad. des Sciences de Pans, ann. 1786 (Paris, 1788), 
pp. 61&-643 and 644-683; Traite des F auctions Elliptiques, Paris, 1825, t. i. 
p. 79 et seq. 

12 



4 NOTE ON LANDEN'S TRANSFORMATION 

The method discovered independently by Gauss and called by him 
the "algorithm of the arithmetico-geometrical mean" originated in 
connection with the evaluation of complete elliptic integrals arising 
from a problem in attractions required in planetary theory*: in a note he 
mentions that the results were obtained by him as part of a more com- 
prehensive theory, independently of the results of Lagrange and Legendre 
with which they are closely connected. Gauss employed a trigono- 
metrical transformation along the following lines. 

Commencing with two numbers (a , /> ), successive numbers (a n , b n ) 
are calculated from the recurrence formulae, 

+! = I ( + &), &n + l^N/(A) .............. (9). 

Writing A n = J(a n 2 cos 2 </>, t + b^ sin 2 < re ), Gauss employs the recurrence 

formula 

<I H tun < w +i = A n +! tan <)> , 

or tan <j> M - tan <K Jfa+iM - %/(* + ^0/N/K + A n ) . . .(10;. 

From the relation 

tan < ft = ( Aj A 2 . . . A n . . . )/( i . . . n . . . ) . tan < 
we have ultimately, since the <'s tend to a limit $ n , 



and hence 

^A = tan- 1 [(A^,, ... A n ..)/Ki - -0 - tan* ] ...(11). 



Extremely convergent series for the elliptic integral of the second 
kind were also obtained by Gauss. The application of this method to 
the numerical calculation of the elliptic integrals was considered by 
Jacobi some years later t. 

In the course of researches carried out between 1797 and 1808, but 
not published until after his death considerably later j, Gauss developed 
the theory of elliptic functions to a remarkable extent and derived a 
large number of formulae which, at a later date, were obtained inde- 
pendently by Abel and Jacobi and which are now associated with the 

* Gauss, "Determinatio Attractions quam in Punctum quodvis Positionis datae 
exerwret Planeta, etc.," Comm. Gott. Soc. Reg. Scient. iv. 1818; Ges. Werke 
(Gottingen, 1866), Bd. in. pp. 331-355. The note mentioned above is dated Feb. 9, 
1818 ; Ge*. Werke, Bd. in. pp. 357-360. 

f Jacobi, Fundamenta Nova, 1829, sections 38 and 52 ; Ges. Werke (Berlin, 1881), 
Bd. i. pp. 154 and 203. Also " Numerische Berechnung der Elliptischen Functionen," 
Crflb; Bd. xxvi. pp. 93-114; Ges. Werke (1861), Bd. i. pp. 345-368. 

J Gauss, Ges. Werke (1866), Bd. in. pp. 361-403 and 433-469. 
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theory of the Theta-functions. In one of these posthumous papers, Gausi 
employs the recurrence formula tan20 n = N /(& w +i/ a n+i)tan< n4 . 1 , which 
is much more convenient for purposes of numerical calculation than (10). 
The last mentioned method is considerably simpler than that of 
Legendre's scales of moduli, especially when a calculating machine is 
available, although in theory these two methods of computation are 
practically identical. 

Although a few formulae based on the scales of arithmetico-geometrical 
means were briefly considered by Jacobi t, no attempt seems to have been 
made by him to apply the results to actual computation. 

As the method of evaluating the elliptic functions by the use of the 
A.G.M. scales is only briefly dealt with in existing text-books t, and does 
not seem to have been developed to the extent it deserves in various 
memoirs which have appeared on the subject in recent years, a tolerably 
full account of the subject is given in the present book, together with 
a brief outline of the main analytical steps. Several new formulae have 
been obtained by the writer which facilitate the application of this 
method of numerical computation and considerably extend its scope. 

From the point of view of pure analysis, the modern Weierstrassiari 
notation presents many advantages in regard to elegance and symmetry, 
a superiority which no longer holds, however, when it is required to 
proceed to numerical evaluation. As the main point of the present 
paper is concerned with improvements in methods of computation, the 
older Jacobian notation has been adhered to as far as possible, especially 
as it expresses most readily the elliptic functions and integrals which 
arise from integrations involving circular and hyperbolic functions. The 
results are set out in such a manner as to provide a fairly complete 
compendium of formulae intended to be of service to the computer in 
reducing to numerical results such expressions involving elliptic functions 
and integrals as may arise in physical and astronomical problems. 

* Gauss, loc. cit. p. 388. 

t Jacobi, Ges. Werke (1881), p. 357. 

$ Cayley, A., Elliptic Functions, 2nd Ed. (G. Bell and Sons, London, 1895), 
pp. 326-338 ; Enneper-Miiller, Elliptische Functionen (2nd ed. 1890), pp. 361-364 ; 
Tannery-Molk, Fonctions Elhptiques, t. iv. 1902, pp. 269-275. According to a foot- 
note in Halphen's Fonctions FMiptiques, t. 11. 1888, p. 310, the subject of the 
arithmetico-geometrical mean was to have received special attention in the author's 
third volume, unfortunately incomplete at the time of his death. 
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SECTION III 

ON THE SCALE OF ARITHMETICO-GEOMETRICAL MEANS 
In forming the scale of arithinetico-geometrical means, we start with 
the positive numbers (a 6 ), of which a is the greater, and form 
successively 



...(12). 

* > 

We obviously have the relations 

The a's and 6's tend to the same limit denoted by 



(14) 



with extraordinary rapidity, even when a and ^ are initially numbers 
of very different magnitudes*. 

We notice in passing that M ( , 60) satisfies a homogeneity relation 

of the form 

M(ai,b ) = M(*a ,d> ) ................... (15), 

being any number. The array of numbers (12) will be referred to as 
the " scale of arithmetico-geometrical means (a , J )" or more briefly 
as the A.O.M. scale (a , 6 ). Also, the limit M(a 09 b ) will be denoted 
by a n , as long as by so doing no ambiguity is involved. 

If we calculate in the same way an array of numbers commencing 
with {a ' = a , A ' = c = V(o a - V)}, we derive what will be called the 
"complementary A.O.M. scale (a '> 6 ')," the symbols being in this case 
denoted by accented letters, and the limit M(aJ, &</) by #'. 

* The usual arithmetical process of extracting square roots is easily adapted to 
most calculating machines of modern type. An account of these is given in 
Modern Instruments and Methods of Calculation, edited by E. M. Horseburgh 
(London, G. Bell and Sons, 1914); also by d'Ocagne, Le Calcul SimplijU (Gauthier- 
Villars, Paris, 1905). In carrying out the process the work is considerably shortened 
by keeping in mind the rule that if the first p digits out of the number n 
required in the square root have been obtained by the usual process, the next 
p-l digits can be obtained by division only, with a possible error of 1 in the last digit 
(Chrystars Algebra, Part i. 5th ed. 1904, p. 210). To illustrate the rapid convergence 
of the A.O.M. scales, Gauss works out an example starting with a =l and b =0'2 : 
<z 6 and 6 6 are represented by the same number to fifteen significant figures. 
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It is obvious that the operations of the A.G.M. scale may be carried 
out backwards : if the entries are denoted by negative suffixes, they 
are connected with those of the complementary scale by the relations 
a. n = 2 n a n ', b- n = 2 n c n ', c. H = 2 n b n ' (16). 

These relations evidently hold good when accented and unaccented 
symbols are interchanged. As a result of these relations it will be 
seen that there is no advantage to be gained from considering any other 
than the A.G.M. scale of positive suffixes (a n , 6 n ) and the corresponding 
complementary scale (a n ' t 6 n '). 



SECTION IV 

LANDEN'S SCALE OP INCREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second 
kinds, modulus k. 

If we write A w = J(a n * cos 2 < 4- b n 2 sin 2 ^), the recurrence formula 

tan (<K+i - <) = (bnlcin) tan <t> n (17) 

may be written in either of the forms 
sin (2< w - <K+i) - (tWiMi+i) sin < n +i, 

cos^-c^O-A^K+i (18). 

From these we derive 

A+i + c +1 cos +! = * and A n+1 -c n+1 cos^ n+1 -a n ^ n /A n ...(19), 
which give on differentiation, 

^ /A = ^^A 1 -i(^ 2 /A 2 ^...^(i)^ n /A w - (20). 

From the relation 

A 2 n+1 + c n +! A n+1 cos <t> n+l = \ A n 2 + \ a n b n 
we obtain, making use of (20), 

^ n d<l> n - & n +id<t> n+l =C n +i COS < B +iC?< n +l ~ n^n^n/A n ...(21). 

Subtracting from each side the identity 

a*d$ n l n - a\ +l dfa+J^ = } c tt 2 ^ n /A w - a n b n d^ n /^ n . . .(22) 
we obtain 

(A n - O^/AOrf*. - ( A ^l ~ a 2 n^l/^n + i)^n + l 

- C n+1 COS ^ + irf^ + i - iCiWn/A* . . .(23). 

Writing down a series of such equations, commencing with n = 0, 
and noting that as n increases, (A n - a n 2 /A n ) tends to zero in the limit, 
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we have, on integrating, 

/*o r* 

I A d< -a 2 1 d</>o/A = c 1 sin< 1 + c 2 sin<2 + + c n sin < w -t- ... 
./o ./o 

- 1 (c 2 + 2d f + 4c 2 2 4- . . . + 2V n 2 + . . .) f^ /A ..... (24). 



./o 

If we now construct the A.G.M. scale (#0 = 1, &o = #') an( ^ calculate 
successive amplitudes commencing with ^0 = ^, we have A = A(<, k). 
As w increases, the a's and //s tend to the same limit a n , while the 
angle < n /2 w tends to a finite limit. We thus obtain from (20) and (24), 

-.P(* *) = (!/) (^2") ................. (25), 

fi(<t>> *) ~ -P 1 (& *) = fa sin <k + c, 2 sin < 2 + . . . + C M sin ^ + . . .) 

-i(co 2 ^2c 1 2 + 4c 2 2 +... + 2 w c n 2 +-..)^ T (*^) ...... ( 26 )- 

If we commence with < = |TT, we have 

<i = T, ^2 - 2?r, </> 3 - 4?r ... <^ tt /2 n = TT, 
so that K-\TT\a n , and 

(^-A t )//r-i(^ + 2f: 1 2 + 4^ 2 +... + 2V n 2 +...) ...... (27). 

The series (24) and (27) converge with extreme rapidity, resembling 
in this respect the (/-series with which they are closely connected*. 
Jacobi's integral Z, (<, k} is defined by the relation 

Z(h *) = J0(* K)-(EIK}F(<}>, k} ......... (28) 

or, if the angle <f> is connected with the argument u by the relation 

sin < = sn (u, X')i 



........... (29). 

Jo 

Prom (26) and (27) we derive the expansion 

Zfa *) = Ci sin ^i + c 2 sin ^> 2 + +c n mi<j> n + ...... (30). 

In terms of Jacobi's Theta-function (u, k), 



. *) = log [(ii, *)/0(0, *)] ...(31) 
where (0, *) = ^WT/ir) ..................... (32). 



* These relations are discussed by Legendre (Trait6 des Fonctions Elliptiques, 
t. in. Deuxieme Supplement, 1828, p. Ill et seq.). On the connection between the 
A.O.M. scale and the modern transformation theory of the Theta- functions, see 
Tannery and Molk, Fonctiom Elliptiques, t. iv. 1902, note 2, pp. 269-273. A develop- 
ment of the theory of the Theta-functions based on Legendre's scales of moduli is 
outlined by Bichelot, Correspondenz mit Herrn Professor Schroter (K6nigsberg, 
1868). See Examples 19-24 of the Appendix. 
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From (20) we have d</A(0, ) = Q) n d< n /A n , so that the series (30) 
enables us to write 

rZ (</>, k) e?</A (<, K) = 2 ( c n/2 w ) I sin < w c? n /A n . 
i JQ 

On integrating and making use of (19) and (18), we obtain 

r0n 

i sin </> n d<l> n /& n = (l/Cn) log (a n -i/& n ~i) = (l/c n ) log sec ( 

Jo 

leading to the expansions* 



or 

log [0 (N, *)/0 (0, *)] - } log K/A ) + \ log sec 

+ |logsec(20 1 -0 2 )4-...+(^+ 1 logsec(2^_ 1 --^ l ) + (34). 

In the series (34) it is more convenient to write 

/A - cos (0! - < )/eos 
for purposes of numerical calculation. 

(ii) Calculation of the Jacobian Functions sn (M, A), en (u, A'), dn (w, A) 
etc. in terms of the argument it. 

When the argument u is given and the modulus k is known, it is 
only necessary to compute the A.G.M. scale ( =1, b v = k') to such a 
value of n that c n is less than the small quantity determining the order 
of accuracy of the calculations. Since u - /' 7 (<A> k), equation (25) enables 
us to calculate in circular measure the angle < from the relation 

* = 2"ai (35). 

From the recurrence formula (18), 

sin (20 W _! - n ) = (c n /a n } sin <, 

we are enabled to calculate successively the angles < n _i ... < 2 , fa and 
finally , in terms of which 
sn (u, k} = sin <f> , cn (u, k) = cos </> , dn (u, k) = A = cos /cos (0, ) 

(36). 

Z(u y k) may then be computed from the series (30), and <>(#, k) from 
formulae (32) and (34). 

* Formulae (30), (33) and (34) are practically the same as those given by Legendre 
(Trait^ des Fonctions Elliptiques, 1. 1. section 90), and are contained implicitly in 
formulae given in Gauss' posthumous papers. Formula (34) is equivalent to that 
given by Jacobi (Crelle, Bd. xxvi. pp. 93-114; Oesammelte Werke (1881), Bd. n. 
p. 357). 
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SECTION V 

THE HYPERBOLIC SCALE OP INCREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second 
kinds, modulus k 1 . 

Writing i<f> n for <f> n throughout the formulae (17) to (24), we obtain, 
instead of (17), 

tanh (<t> n+l - <f>n) = (b n /an) tanh fa ............... (37), 

which may also be written 

sinh (2 fa - < n +0 = (cn+iA*n+i) sinh ^ +ll cosh (2 fa - <f> n+ J = A n+1 /a TO+1 

...... (38). 

If we denote A n = >J(an cosh 2 fa - 6 M 2 sinh 2 fa\ it is easily proved, as 

before, that 

^ /A = |^ 1 /A 1 = i^/A 2 =... = Q)^ B /A B = ...... (39) 

from which it follows that as n increases, we have in the limit, 



If we now write 

sinh fa = tan < ........................ (40) 

it is easily proved that 

A = A ( #)/cos fa dfa/*o - dfa* (<#>, k'} ......... (41), 

so that 

u = F(+ 9 k') -(1 fan) (<f> n /2 n ) ............... (42). 

Formula (24) becomes, 

/^o r0o 

I A rf^> - 2 I ^ /A = (cj sinh fa + c 2 sinh fa + . . . + c n sinh <#>+) 
Jo Jo 

- J (c 2 + 2r I 2 + 4C, 2 + . .. + 2"c !1 +...)[ ^/Ao -(43). 

Jo 
It is easily shown that 



2 (u, k'} = ( "dn 2 ( 
Jo 

or in terms of the Z- function of imaginary argument, 

I \ Q d<t> = (l/i)Z(iu,fy + (E/ir)u. 
Jo 

Substituting in (43), and making use of (27), we derive the series, 
corresponding to that of (30), 

(\li}Z(iu, k) = Cj sinh fa + c 2 sinh fa + ... + c n sinh fa + ... .. .(44). 
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The following formula, corresponding to (34), easily follows by 
integration, 

log [0 (in, *)/0 (0, *)] - i log K/A ) + ( J) log sech (2^.! - *.) 

...... (45) 

in which we may conveniently write # /A = cosh (<i < )/cosh < . 
From (42) we easily deduce the relations, 
(I ft) sn (iu, k) = sn (u, ')/cn (w, k f ) = sinh < 
en (t, A) = 1/cn (M, K) = cosh < 

dn (iu y k) = dn (w, #')/cn (u, k') = A = cosh < /cosh (^ - < ) 

...... (46). 

Making use of Jacobi's well-known relation in the form 

Z(u 9 tT) = sn (w, k'} dn (w, i^/cn (^, ^) - iru/(2JEZ" ) -f tZ (tw, 4) 

...... (47) 

we have from (28) and (46) 



= sinh < /cosh (& - <#> ) - * 
Noting that Legendre's relation gives 

EIK + EIK'-l=\*l(KK i ) ............ (48) 

we obtain from (27) 

E(h k'} = sinh </>o/cosh (*, - <#> ) + (1 - E/K) F(4>, K) + 1 ^(111, *), 
or, finally, making use of (27), (42) and (44), 
E(<t>, *') - sinh ^/cosh (^, - </> ) 

+ \(ef + 2 Cl 2 + ... + 2*^+ ...) (IK) (*/2 rt ) 

- (d sinh <#>! + c 2 sinh <^ 2 + + c n sinh <f> n + . . .) .(49). 

(ii) Calculation of the Jacobian Functions sn (w, k'}, cn (w, 4 ; ), 
dn (u y k'\ etc. in terms of the argument u. 

Compute the A.G.M. scale ( = 1, b = lf) and calculate <f> n from the 
relation </> w = 2 w a w w, w being such that c n is less than the small quantity 
determining the order of accuracy of the calculations. The hyperbolic 
recurrence formula 

sinh (20,*-! - <) = (Cn/aO sinh </ 

enables us to calculate successively < n _i ... < 2 , <^i and finally <<>, in 
terms of which we have from (46) 

sn (u, If) = tanh < , cn (u, k') = 1/cosh < , dn (w, A') = 1/cosh (^ - </> ) 

...... (50). 

may then be calculated from (47) and (44). 
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(iii) Calculation of the Jacobian Functions of imaginary argument, 
sn (iu y ), en (, A), dn (fti, i), 0fo. 

The procedure outlined above, together with equations (46), obviously 
enables us to calculate the Jacobian functions of imaginary argument 
and modulus , while formulae (44) and (45) serve for the functions 
Z(iu, k) and 0(i, k). 



SECTION VI 

GAUSS' SCALE OF INCREASING AMPLITUDES 
(i) Circular Hcale. 

In connection with the A.O.M. scale, Gauss made use of the recurrence 
formula* 

tan Xn+i = N/(n+ Ai+i) tan 2xn ............... (51). 

It is easily proved that the above formula is connected with Landen's 
scale (17) by the relation 

tan Xn = ,J(b n la n } tan < n ..................... (52) 

or 

*.+i-* = 2x-i ....................... (53). 

In the limit, as n tends to infinity, (b n fa n ) approaches unity, so that 
Xn and < both tend to the same limit. It thus follows from (25) that 
if we form the A.G.M. scale (a =- 1, & ( ,^X;'), and starting with Xo deter- 
mined by tanxo - \/(^>A*) ton <> compute successively Xi Xs> X 
we have 



From (52) and (53) it is a simple matter to derive successively 
<, , < 2 . . . < n . . . for use in the formulae of Section IV. When these angles 
are required, it is, however, simpler to make direct use of Landen's 
recurrence formula (17) or (18). 

(ii) Hyperbolic Scale. 

It is evident that formulae (51) and (52) may be written in the 
hyperbolic form 

tanh Xn+i-N/(^W ft n+i) tanh 2x* ............... (55), 

tanh X n = >/(&AO tanh < n .................. (56). 



* Gauss, Ges. Werke (1866), Bd. in. p. 388. To compare Gauss' results with 
those of the present paper we note that W =2 II F H , x n =2 n ^ n " i " 1 - 1* t^ en appears 
(p. 392) that formula (26) of the present book was also obtained by Gauss, inde- 
pendently of Legendre (foot-note on p. 9). 
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As n increases indefinitely, Xn and < tend to the same limit, so that 
if we write 

tan< = sinh0 or sin < = tanh < ............... (57) 

it follows from (42) that 

F(<t>,^ = (l/a n )(t,J2"^(l/a n KXn/2 n } .......... (58). 

Noting that, as before, 

<*+!- 4>n = 2x-i ................. (59) 

it is a simple matter to derive successively 0,, </ 2 , ... <f> n ... for use in 
the formulae of Section V. In such cases, however, it is simpler to 
make direct use of the hyperbolic recurrence formulae (37) or (38). 
(iii) Calculation of K\ E' and the nomc q~-e -"^'/A" 
If we write < = JTT, it follows from (57) that tanh <,, 1, and from 
(56) that tanh x = V( AO- ^om (55) we deduce that 

tanh 2xo =&i/i, tanh 2xi = /> 2 /2, tanh 2x =- 'VnA'n+i, 
and thus 2x - log {(a n+l + * f ,)/(+i - '> w +i)} - 1 log (tf^AW;) 
or (x/2 w )-(i)^ 2 logK, 2 A',, +a ). 

It follows from (58) that as n is indefinitely increased, (i)" log(ti n /c n ) 
tends to a finite limit given by 

<i K F(\*,k') = ^K'IK- J log (1 A/) --.Q)" log (*./<?,,) (60). 
Making use of the relation (13), dc n - J^ 2 M-1 , we have 

log (/<) - 2 log (2a n /a n ^) + 2 log (WH.,/^-,), 
which gives by successive application, 



2 log (2^/^-0 + 4 log (2u n _,/W w _ 2 ) + . + 2 n log (2a,/f/ ) + 2 n log (fl ffl /r ). 
Inverting the order of the series, we obtain 



+ log (0,K) + i log (a 2 /a,) + .. f (J)"- 1 log (/_,). 
As n is increased indefinitely we obtain the series 

rt n /T - i log (4^/0- i (J)" log^/^/^O ...... (61). 

Making use of Legendre'.s relation (48) in the form 



and using (27), we obtain the following useful formulae, 

) (^o 2 + 2c," 4- 4r, 2 + + 2cv + - - ) 
x {log (4/1,/r,) - log (aja >) - i log (^,/a,) - }-* ...... (62), 

= Kl/)(l-ko 2 -^-2r 2 2 -.. -2"-V--0 

x {log (4j/Ci) - log (tf,/^) - J log (tf a/tf,) --..}- ...... (63) 
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for the calculation of the complete elliptic integrals of complementary 
modulus in terms of the A.G.M. scale (a = 1, i = *') 

The nome q is denned by the relation log (l/q) = 7rK'/lC, from which 
it follows that 

log q = 2 ( i)- 1 log KK+0 - log (40,/d) (64) 

which is more convenient for numerical computation than (60). It may 
be noted that the entire series of calculations may be carried out with 
ordinary logarithms to the base 10, enabling q and powers of q to be 
evaluated directly*. 

The complementary nome q' - e~"" K l K ' is obviously related to q by the 

formula 

log q . log q = 7T 2 . 

SECTION VII 

LANDEN'S SCALE OP DECREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second kinds, 
modulus k'. 

If we denote 

^---J^cos^n-t-^sin 2 ^) (65), 

the recurrence formula 

tan ($ n - I/VH) = (c+i/a n4 .i) tan I/VM (66) 

may be written in either of the forms 

sin (2^+i - i/O - (&/) sin ^ n , cos (2^ H+ i - ^n) = A ;l K . . .(67). 
From these we deduce 

A n + ft ll cos^ H -=2A n+1> A n -6 n cos^ n --- 2rt n+1 6' n+1 /A ro+1 ...(68), 
which give by differentiation 

ify /A = rfft/A, = <fy. 2 /A 2 = . . . , <fy n /A n - (69). 

If we construct the A.G.M. scale (a - 1, b^ = k') we notice that as n 
increases, \f/ n tends rapidly to a finite limit. If we start with ^ ^</s 
A = A (^, k'\ we have, since r n ultimately tends to zero, 

* The limit derived in (60) does not appear to be generally known, although an 
equivalent formula is given by Legendre (Fonctions Elhptiques, t. i. pp. 94, 101). 
Although formula (GO) does not appear to have been explicitly stated by Gauss, 
the scries (61) is found among a number of formulae given in his posthumous 
papers (Qauss, Gen. Werke, Bd. in. p. 377). Neither (60) nor (61) appears to have 
been noticed by Jacob! or by subsequent writers on elliptic functions. 
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^ 
') = (1/O / n ^/cos^ B 

J 

a useful formula due to Legendre*. 
From (68) we easily obtain the relation 

A n 2 + b n A w cos ^ - 2A 
which gives 



Writing down the system of equations for n - 0, 1, 2, ... 71, multiplying 
successively by 1, 2, 4, ... 2", adding and integrating, we obtain the 
result 

+i rfo 

A n+1 ^ n+1 = i(c 2 4.2^ 2 -i-4c 2 2 4-...+2 w O / <%/A 
o /o 

- (6 sin ^o + 26, sin ^i + ... +2 n 6 n sin^ IB ) ...... (71). 

Adding to each side the identity 



,+! cos </Oi+ 1*^11+1 = 6 n +! sin i/^ n+1 + (2 W " M 1) 6 n+1 sin ^r n+1 
and re-arranging the terms, (71) may be written 

[fa f^n+l 

Jo yo 

/"\^0 

= -J- (c + 2^ + . . . + 2 n c n 2 ) I d\l/J& + 6 rt +i sin \I/ n+l + (b l sin \!/ l - 6 () sin \l/^\ 
JQ 

+ 3 (6 2 sin ^ 2 - bi sin ^) + . . . + (2 n+l 1) (6 rt+1 sin ^ t+1 6 W sin \f/ n ) 

(72). 

It is not difficult to establish the identity 

"(73). 

As n increases indefinitely, On+i and 6 n+l tend to the same limit a n , 
while c n+1 tends to zero. It follows that the second integral on the left- 
hand side of (72) tends to zero, so that we finally obtain, making use 
of (70), 
E (i/s k r ) = a n sin \l/ n + 1 (c 2 + 2c* + . . . * 2" r rt 2 + . . . ) 

x (l/a^ log tan (Jir + \\l/ n ) + 2 (^ , c n ) (74), 

2 

where for the sake of brevity we denote 

2 (\l/ n , c w ) 2r -x " - 2 - . ] . + 6c 3 x tj 2 N 4- . . . 

2 COS (2^j ^ a ) COS (2^4 \f/ A ) 



* Legendre, Fonctiom Elliptiques, t. i. chap. xix. 
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a useful fonnula for the numerical calculation of the second elliptic 
integral*. 

Making use of (27), (28), (48) and (70), we obtain from (74) the 
following interesting expressions 

Zty, A') + \*J?(*, k'}l(KK'} = a n sin ^ + 2 ty n , c) ...(76), 

2 

or, in terms of u, by (70) we have 

sini/^ rt = tanli (a n u\ } 

and Z (u, k'} + a n uj K' = a n tanh (a n u) + 2 (i// n , C M ) f ' '( 77 )' 

Jacobi's relation (47), taken with equation (76), gives 

(I/O Z (iu, /') - tan ^ cos (2fc - ^ () ) - a n sin ^ - 2 tyn, f (78). 
The addition theorem 



v ...... (79), 

with iu for u and v = (K+i/T) gives two interesting formulae 
iZ(iu + ^T+ lA", /) = (!+ sin M ) 




If we notice that 

/"^ H 

I sin ^//^M/AH = (!/&) log (rt,n.i/A+i) ^ 
yo 

the result of integrating (72) written in the form 

00 

Z(u, k') + a u u//C = <9 H sin \l/ n + 2 (2 n+1 -l) (b M sin ^ w+l - b n sin ^ n ) 
gives rise to a formula for the numerical calculation of (u, &'), 



in which, if desirable, we may write sec ^ n ^ cosh (<?/). 

(ii) Calculation of the Jacobin n Functions sn (?/, /'), cn(w, X y ), dn (w, /:'), 
^c % . /w fc/v/*tf r>/' #//<' argument it. 

Form the A.U.M. scale (<i,,-l, (>o=~-k'} and determine ^,, from 
sin ^,, tanh (<///), w being such that c n is less than the small quantity 
determining the order of accuracy of the calculations. 

* The series (71) was obtained by Legendre (Fonctions Elltptiques, 1. 1. chap, xxi.) 
and left by him in this non-com ergent foim. As far as the writer is aware, the 
convergent series (75) is new. 
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From the recurrence formula (66), tan (i/r w - 1 -^ w ) = 
compute ftt-1 -"ft, ^i, and finally ft, in terms of which 
sn (u, V) = sin ft, en (u, #) = cos ft, dn (u, ') = A = cos (2ft - ft) 



may then be calculated from (77) and (tf, A') from (81). 



SECTION VIII 

THE HYPERBOLIC SCALE OF DECREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second kinds, 
modulus k. 

If we write t^ for ^ throughout the formulae (65) to (69), (65) becomes 
A n =V(n a cosh 2 ft>-c tt a sinh 3 ftO ............... (83), 

while the recurrence formula (66) may be written 

tanh (ft> - ft, +1 ) = (c n+1 /a n+1 ) tanh ft, +1 ......... (84), 

or, more conveniently, 

sinh (2ft+! - iAn) = (&/) sinh ^ n , cosh (2^ n+1 - iAn) = A n /a n . . .(85). 
As in Section VII (i), it is readily proved that 

^ /A = ^ 1 /A 1 = ^/^a-...=^n/A n = ......... (86). 

If we construct the A.G.M. scale (a = 1, b = fc') and write 

tan^ = sinh ^ ........................ (87), 

it is readily proved that 

A = A (^, ^)/cos $ and rfi/r /A = d/r/A (ft *) ...... (88). 

Thus, writing w = J^(ft ^), we have 



sn (u, K) = sin ^ = tanh ft, en (u, K) = coa\l/ = sech ftl . , 
dn (M, k) = A (ft A:) = A /cosh ft = cosh (2ft - ft>)/cosh ftj ' ' ^ '' 
As n increases, the ^s tend to a limit ft, while c n tends to zero. 
Thus the result of integrating (86) gives 



u = JF(ft *) = (l/a n ) rfft/cosh^ n = (l/a n )tan- 1 (sinhft) ...(90). 
If we notice that 



f dn a (u, *) <fo/cn" (M, *) = 
JQ 
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we have from (29) 



= (1/0 Z(iu, *') + (E'IK"> u. 

It is easily seen that equation (72) holds good when i$ n is written 
for \l/ n . Hence if we write 

r r M 

i/' 8 -^j) 



+ 2 C2 - 1 ) C *- + 

+ 2(2 lK + , cogh(2|/ , B+s _ |/W2) + 

and make use of (47) and (48), we obtain after some reductions 
Z (fr k) = tanh ^ cosh (2fc - <Ao) - <* sinh ^ n - 2' (^ n , c tt ) . . .(92). 

The relations (27) and (28) enable us to calculate the second elliptic 
integral Ety> k). 

(ii) Calculation of the Jacobian Functions sn (M, A), en (u, k), dn (u, k), 
etc. in terms of the argument u. 

Compute the A.G.M. scale a Q = l, b Q = k' and determine ^from (90), 

which may be written 

sinh \l/ n = tan (a n u) ..................... (93), 

n being such that c n is less than the small quantity determining the 
order of accuracy of the calculations. 

We then calculate the amplitudes ^ n ^ ... ifo, <Ai, <Ao from the recur- 
rence formula 

tanh (i/w - i/O - (c n /a n ) tanh ^ n . 

Equation (89) then gives 

sn(w,A;) = tanh^ , cn (u, A) = sech M 
dn(w,^) = cosh(2^-^ )/cosh^ J 

Equations (91) and (92) then enable us to calculate Z(u, k). To 
compute (/, k\ we note that (92) may also be written 
Z(u, k) = sn (M, *) dn (u, A)/cn (u, k) - a n sinh $ n 

sinh ^ w ^ - b n sinh fc). 







Integrating this expression with respect to u and noticing that 

^L 

\ W sinh tAn^n/An = (l/*) log cosh (2^ n+a - ^ +1 ) f 
JQ 

we obtain the following series for (w, *), 



...... (95). 
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SECTION IX 

ON THE NUMERICAL COMPUTATION OF THE 
THIRD ELLIPTIC INTEGRAL 

As defined by Legendre, the third elliptic integral may be written 

............ (96). 



The parameter n is taken as real, and may take any value, positive 
or negative. 

Writing u = JP(<, A), n = - 2 sn a (a, A), Jacobi considers the integral 

. N ,, nS)(a--tt)~| ,^. 
wZ(a) + Hog ) -- : ...(97). 

\ / - & 



Legendre's integral is connected with Jacobi's by the relation 

n s (n, k, <) - u + sn a/(cn a dn a) . IT (w, a) ...... (98), 

and in this book is differentiated from it by the suffix 8 . 

In order to include all real values of n in (96), a is not restricted to 
be real, and the treatment is usually divided into four cases : 

Case I. n negative between and - K 2 ) rr . ,. 
n TT 4.- u x * j c Hyperbolic Cases. 

Case 11. n negative between 1 and <x> J ^ 

Case III. n negative between -#* and - 1) ~ , 
r* TTT "j.- v. j. /A j r Circular Cases. 

Case IV. w positive between and + oo J 

CASE I. ft negative, between anc? /fc 2 . Hyperbolic case. 
We write w = - ^ sn 2 (a, k) = - ^ sin 2 0, where < a < IT . . .(99). 
Prom (97) and (98) we have 



*r / 7 ^\ T, 
H 8 (n, *,<#>)= w 1 + 
8V ' '*' L 



/ 
en (a, 



____ 

Cn (a, *)"dil (a, *) (a + U, k) 

Compute the A.G.M. scale (a = 1> b = k'\ and starting with < = ^, make 
use of the recurrence formula (17), tan (< n +i -#*) = (6 n /n) tan <^ n . We 
then have sn (a,)=sin , en (a, #) = cos <^ , dn(a, K) = cos /cos(^i~^ ), 
and <ZT(a, ^) c i si n 0i + ^2 sin <^ 2 + + CH sin <^ n + . . . , while a n a = 
We thus have 



22 
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The calculation of u by the recurrence formula (17), commencing with 
</>o = fa presents no difficulty. The most convenient method of computing 
the last term of (101) is to make use of Jacobi's series 

(x) = 1 - 2q cos 2a n as + 20* cos 4a n # - 2q 9 cos 6a w # + ...... (102) 

in which we have written a n =^/K. The nome q is most conveniently 

oo 

calculated from (64), log q = 2 ($) n ~ l log (a n /a n+l ) - log (4^/cO. 

If we write < = |TT, the last term of (101) vanishes, and we obtain for 
the complete elliptic integral of the third kind the expression 

H 8 (n, k, JTT) = K[l + tan fa sec fa cos (fa - < ) fci sin fa + c 2 sin < 2 + 

... + c n sinc/ n +...}] ...... (103). 

We may also make use of the hyperbolic scale of decreasing amplitudes 
and the formulae of 8. The final results are not, however, so con- 
venient for computation as (101) and (103). 

CASE II. n negative, between - 1 and - oo . Hyperbolic case. 
We write 
n = - If sn 8 (a + $ JT, *) = - 1/sn 8 (a, k) = - I/sin 2 0, where < a < K 

...... (104). 

Equation (98) becomes 

du 



n,(fc, *, ^) = sn s o f t du a =M -- . 8I L" n(tt,oH-;jr ; ) ...(105). 
8V f '^ 7o sn 2 a-sn a M cnadna v ' ' ^ J 

We have 

^(a + eJr / ) = ^(a)-t-cnadna/sna-^W/jr ...... (106) 

and %(z^iK'^iBH(z) ..................... (107), 

where log B = Jir JT'/^- Jirt 0/^- 

^T() is Jacobi's Theta- function given by the expansion* 

H (z) = 2#* (sin a n z-<f sin 3az + g 6 sin 5a tt - g 18 sin 7a n -*-) 

...... (108). 

Writing successively (a - u) and (a -f u) for z in (107), and a 4- t'JT for a 
in (97), we obtain 

n (w, a -f tTT) = wZ (a, *) + W en a dn a/sn a 4- J log [^T(a - u)/H(a + tt)], 

and hence from (105), u < a, 

H 3 (n, *, *) = - sn a/(cn a dn a) [uZ(a, k) + l log {H(a-u)/H(a + 1)}] 

...... (109). 

* Jaoobi, Fundani/nta Nova, 62, 63 ; Ots. Werke, 1881, Bd. n. pp. 229-231 ; 
Cayley, Elliptic Functions, 2nd Edition (1895), Chap. vi. 
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For numerical calculation, form the A.G.M. scale (o= 1, <> = # 
starting with ^0 = ^, make use of the recurrence formula (17), 

tan (< n+1 - <) = (6 n /a n ) tan < n . 
We then have a n a = < n /2 n , and as in Case I, we derive finally 



= - u tan < sec # cos (& - < ) fci sin <^ + c a sin < a + . . . + c n sin < w + . . . } 
- tan <t> sec < cos (& - < ) log [ff(a - t*)/^"(a + a)] ...... (1 10). 

Starting with < = and making use of the same recurrence formula as 
above, the computation of u presents no difficulty. As before, the 
calculation of the last term of (110) is best achieved by computing the 
nome q by (64) and making use of (108). 

If we write < = JTT, the last term of (110) vanishes, and we obtain for 
the complete elliptic integral of the third kind the expression 



= - .fiT tan < sec < cos (<& - < ) fa sin fa + c a sin ^ + ... +c n sin<+ ...} 

...... cm). 

When >a, it will be noticed that the integral (105) has a singularity at 
u=a, where the integrand becomes infinite. Atthe same time the logarithm 
of a negative quantity appears in (110), since H(a u)=: H(u a). 
It is necessary in this instance to define the path of integration, re- 
garding (105) as an algebraic integral obtained by writing z = sn u. Near 
the point z = sn a this integral is approximately 

H 8 = sn 2 a/(cn a dn a) . / <fe/(sn 2 a - z 2 ), 

which integrated round a small semi-circle below the real axis gives rise 
to the imaginary term - %vi sn a/(cn a dn a). If we define the path of 
integration of the algebraic form of the integral as consisting of this 
small semi-circle and the remaining portions of the real axis, the above 
imaginary term cancels the term - sn a/(cn a dn a) . log ( 1) appearing 
on the right-hand side of (110) when we write H(o.-u) = -H(u-a). 
With this meaning attached to the path of integration, we may write 
the incomplete integral in the form, u > o, 

n, (, *,<) = - sn a/(cn a . dn a) [uZ (a, k) + | log {H (U - a)/#(w + a)}) 

...... (112). 

Numerical computation proceeds in the manner already described. 
CASE III. n negative, between - ft and - 1. Circular case. 

First Method. 
We write 

==-A 2 8n a (ea + jr,*) = -* a /dn a (a,A') = -.^/(l-*' a sin^) ...(113), 
whence sin = sn (a, #), < a < K\ 
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Writing ia + K for a in (98), we find 

H 8 (n, k,<fi = u + i dn (a, k")/[k'* sn (a, #) en (a, #)] . H (u, ia + JT) 

...... (114), 

while (97) gives 



Making use of the addition formula (79) for Z(u), we find after some 
reductions 

iZ(ia + K) = iZ (ia, k} + If sn (a, yfc')/[cn (a, A') dn (a, A')] . . .(116), 
or, from (47), in terms of a real argument 



...... (117). 

Introducing Jacobi's Theta-f unction j (z) defined by the relations 

a(* + JO=a,() ..................... (118), 

i ()=! + 20 cos 2a w 5J 4- 20* cos 4a + 2^ cos 6a n z + 
we may write (114) in either of the forms 



_ 

an (a, A') : en (a, A') 



For numerical computation we may make use of two sets of formulae, 
according as we employ circular or hyperbolic functions. 

(i) Use of circular recurrence formulae. 

Form the A.G.M. scale ( = 1, fr = *') Commencing with < = & make 
use of the recurrence formula (17), tan (</> n+1 - <) = (/) tan < n , and 
thus derive a n u - <^/2 n = <. 

Commencing with ^ = ^i employ the recurrence formula (67), 
sin(2^ 1 -</r n )^(^ n /a n )sm^ n , giving a n a = logtan(j7r + i/r tt ). The 
application of (47) and (76) to both (120) and (121) gives us the formula 



...... (122) 

for the numerical calculation of the incomplete integral, and the formula 
H, (n, *, |ir) = 2 GST/A") cos (2^ - ^)/sin 2^ . [. sin f B + 5 f> )] 



NUMERICAL COMPUTATION OF THIRD INTEGRAL 23 

The peculiar difficulty arising in the numerical computation of the 
circular cases of the incomplete elliptic integral of the third kind lies in 
the evaluation of the T heta-functions of complex argument. Writing 
2 = tt+fa in Jacobi's expansion (119), the series may be resolved into 
real and imaginary parts*. If we write 

8in2a u . htan+ ty 4 sin 40 n u 



_ ______ 

~ 1 + 2q cos 2a n u . cosh 2a n a + 2^* cos a n u . cosh 4a w a + . . . 

...... (124), 

we easily find log[,(w-eo)/0, ( + i'a)] = 2t, from which it follows 
that (122) may be written 

H 8 (, #, </>) = (2/#) cos (2fc - V )/sin 2^ [M {a K sin ^ + 2 (A)} - *] 

2 - ..... (125). 

Powers of the nome q arising in (124) are most conveniently calculated 
from (64), 

log q = 2 Q) 7 *" 1 log (a n /a n+l ) - log (4 I /c? l ). 

The use of hyperbolic functions in (124) may be avoided by noting 
that if we write 



2 cosh 2a n a = (+ Q~\ 2 cosh 4a n a - Q 4 + 
2 sinh 2a n a = Q 2 - Q- 2 j f 2 sinh 4a n a - (^ - 

from which it follows that (124) may be written 



...... (124 to). 

(ii) Z7s of hyperbolic recurrence formulae. 

Write = -A 2 sn 8 (ia4-JST,A;) = ~A 2 /(l-*' 2 tanh 2 ^o / ) ...... (126), 

whence sn (a, A') = tanh <', < o < ^T'. 

From the hyperbolic recurrence formula (37) 

tanh (<// ftfl - */) = (6K) tanh < n ' 
we have 

aa = *n72 tt = ^', sn (a, #') = tanh <^ ', 

en (a, A') = sech ^ ', dn (a, A') = sech (&' - ^')- 

* This procedure is doe to Legendre, Fonctiont Elliptiquc8 t Deuxieme Supplement 
(1828), iz. p. 146 ft eq. 
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Thus from (44) and (120) we obtain 
* C08h +1 



To calculate u we make use of the recurrence formula (17) 

tan (<+! - </> n ) = (6,,/flhi) tan <, giving aw = <k/2 w = $. 
Equation (124) then becomes 

2# sin 2$ sinh 2$' + 2/y 4 sin 4$ sinh 4<E' + 2q* sin 6^ sinh 64>' + . . . 
"" l + 2# cos 2$ cosh 2$' + 2</*cos 4<^ cosh 4^' +2^ 9 cos 6* cosh 6^'+ ... 

...... (128). 

Equation (127) then gives 



**S \. ftf i "'I YV / /a A 1 i/ 1 

8v > -r/ Ar a tanh0 cosl 
and 



TT / 7 i \ -^ cosh <</ f s i nn ^o' S% u JL '1 

n, (n, , -^ TT) = TT; r~iTT> "TV"L / JL ^ wT^A^"^^ 1 " 1 " 1 ^ 
v >2 x A: tanh cosh (^/- ) Lcosh (0! - <^o ) i J 



Second Method. 
If we write = K' - a, (1 13) takes the form 



...... (130). 

We easily find 

f , , N dn (/?, ') 

- 



where J?i () is Jacobi's Theta-function defined by the series 

Hi (z) = 2gi (cos a % z + ^ cos 3a w s + g 8 cos 5a n + q 13 cos 7a tt -...) 

...... (133> 
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Writing V^iloelHifr + iftlH^u-ip)] ......... (134), 

we find 

tan *' = s ^ na ~^ s ^ n h^n^+y a sin3 n Msinh3a n /g-t-<y g sin5a n 8inh5ff n ff+... 
cosa w wcosha n )8+^ a cos 3a n ttcosh3 n ^+^cos5aMCOsh 50^7^+ 

...... (135). 

The numerical computation may be carried out by two methods. 

(i) Use of circular recurrence formulae. 

Starting with ft = # in (130), make use of the recurrence formula (67), 
sin (2ft 4l - $ n) = (in/On) sin ft. Prom (80) we have 

iZ (K+ iK' - iff) = a* (1 - sin ft) - 2 ty), 
while 

sn (ft ') = sin ft, en (ft i') = cos ft , dn (ft 4') = cos (2ft - ft). 
Thus (131) becomes 

...(136). 



As before, u is most conveniently calculated from the recurrence 
formula tan (<t> n+l - < n ) = (b n ja n ) tan <fc t , giving a a u = </ n /2 n = <fc. 

The use of hyperbolic functions in (135) may be avoided by writing 
log Q' = a n /3 = log tan (Jir + Jft), so that 

2coBh B jB=Q'+Q'"M 2 cosh to.j8 = Q +</-! . 

2 sinh a n p = q- Q'- 1 / ' 2 sinh 3 n ft = Q' 8 - (/-'I ' ' ' ' 

Thus (135) may be written 

tan *' = 



. 

Q'" 1 ) cos ^> + ^ (Q'* + Q'-*~) cos 3* + ^ (V s + ^'- 5 

...... (138). 

Writing u- JTin (136), V = TT, and we have, since ^ir^ 



(ii) Use of hyperbolic recurrence formulae. 

Write sin 6' = sn (ft k') = tanh <<>', and (130) becomes 

Make use of the hyperbolic recurrence formula (37) 

tanh (</>'+! - &/) = ( 
and we have 

a n ft = <fr' = < tt '/2*, sn (ft *') = 1 
en (ft K) = sech <</, dn (ft #) = sech (ft' - ^ ')- 
Also - iZ(ift, k) = G! sinh ft' + c 8 sinh & + . . . + c n sinh <^ n ' + . . ., 
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so that (132) gives 
n. (,**) 

_l __ coshc ' 



x f- utf sinh <' cosh (&'- <<,') + u 2 c n sinh^ n ' + *'] . . .(141). 

As before, we calculate u from the recurrence formula 

tan (<f> n+l - <) = (6K) tan <, 
giving a n w = n /2 w = $. Equation (135) gives for *' 

t V - 8mn 3*' g i n *** + ff 2 smn 3*' sin 3<S + <f sinh 5<fr'sin 53> +^ 
cosh $'cos $ + (fcosh 3$' cos 3<3> -f /cosh 5$' cos 5^> + . . . 

...... (142). 

Writing u = K in (141), we have *' = JTT, and since J?r = 



...... (143). 

CASE IV. n positive, between and oo . Circular case. 

First Method. 

Write rc = - A 2 an 2 (to, 4) = JP sn 2 (a, A')/cn 2 ( flj jf) = >p tan 2 tf . . . (144), 
giving sin0 = sri(a,#), 

Thus 

Sn 



...... (145). 

(2;) is Jacobi's Theta-function defined by the series 

(z) = 1 - 2q cos 2a n s + 2^ cos a n z - 2^ cos Ga n z + . . . (146). 
Writing * = f log [0 (w - ta)/0 (u + ^a)] ............ (147), 

we find 

tan * = 



-- _ 
1 - 2q cos 2a n M cosh 2a n a + 2^ cos 4a u cosh 4aa - ... 

...... (148). 

The numerical computation may be carried out by two methods. 

(i) Use of circular recurrence formulae. 
Make use of the recurrence formula si 
commencing with $ Q = 6 in (144). We then have 

sn (a, A') = sin ^ , cn (a, K) = cos ^ , dn (a, A') = cos (2^ - ^ ), 
and by (78) 

iZ(ia, k) = a n sin * + 2 ty) - sn (a, *') dn (a, #)/cn (a, tf) ...(149). 

2 
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Hence (145) may be written 



...... (150). 

As before, u is calculated from the recurrence formula 

tan (<t> n+l - <f> n ) = (bjan) tan < n , 

giving a n w = </> tt /2 n = 4>, while aa = logtan(j7r + i/> n ). We have 
2cosh2a tt a = Q a +Q-M 2 cosh 4a n a = Q* + Q~M , 

28inh2a n a=^-Q-*/' 2sinh4a n a= <*--< ' elc ....... ^ 

where # 

so that (148) becomes 



tan *= 



- 2 )cos 2* i J* -i- Q-*)cos 4* - ^ 9 (0 s + Q L6 )cos6*- . . . 

w ...... (152). 

Writing u = JSTin (150), * = 0, and 



...... (153). 

(ii) f/5^ of hyperbolic recurrence formulae. 

Write sin^ = sn(a J A:')-tanli<^ / ................ (154), 

so that w = - 2 sn 2 (to, *) = If tan 2 tf = If sinh 2 ^'. 

Starting with < ' given by (154), make use of the recurrence formula 

tanh^'n^-^nO-^n/^tanh^' ............ (155). 

We then have 

a<* = 4n'/2 n = *', sn (a, if) = tanh ^ 
en (a, #) = sech ^ ', dn (a, K) = sec 

CO 

and - iZ (eo, A) = X c n sinh ^ n '. 

Equation (145) may then be written 



As before, w is calculated from the recurrence formula 

tan (<+! - < n ) = (ftn/Om) tan < n , 
giving a w w = < n /2 n = <l>. Equation (148) then gives for * 

2g sinh 24> ; sin 2O - 2^ sinh 4<t>' sin 4^ + 2/ sinh 6*^smj6* - . 
"" 1 - 2g cosh^i' cos 2^+ 2^* cosh 4^'cos 4^>- 2^ 9 cosh 6*' cos 6* + .. 

...... (157). 



28 NUMERICAL COMPUTATION OF THIRD INTEGRAL 

Writing u = K in (156), * = 0, and 

* (ll| 4 iW .jr[^ 

Second Method. 
If we write p = K'a in (145), we have 

n - - A 8 sn a (iK* - i/3) = en 2 (5, #')/sn 9 (0, #) = cot 2 . . .(159), 
giving sn (,#) = sin 0, 0</?<JT. 

We easily find 

,<A) 

sn (/?,') en GM' 



where -T(^) is Jacobi's Theta-function defined by 

/T(5?) = 2g* (sin a n z-(f sin 3a n # 4- ^ sin 5a n 5? - 

...... (161). 

If we write 

* = iflog[r(ii + 8)/J5r(-ii + tj8)] ......... (162), 

we find 

, , _ co8a n ^sinha n ff - <? cos_3a n ti sinh 3aJ3 + ^cos 5a n u sinh 6a n P- . . . 
~ sin a n u cosh a n ft - ^ a sin 3a n u coshSa^-f g*sin 5a n wcosh5a)8- ... 

...... (163). 

The numerical computation can be carried out by two methods. 

(i) Use qf circular recurrence formulae. 

Commencing with i^ = in ( 1 59), make use of therecurrenceformula (67) 

sin (2i/wi - ^ w ) = ( 
From (78) we have 



while sn()3,A') = sin^ , cn(ftA;') = co8iAo, dn (ft A') = cos 
Equation (160) then gives 



08(2 ^ 

...... (165). 

As before, u is calculated from the recurrence formula 

tan (0 n41 - < n ) = 
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giving aM = < n /2" = *, while o)8 = log tan ({* + ^,). Thus 
2 cosh , = # + # 



= "-- ' etc ....... 

where $ = tan ( J ?r + 

Thus (163) becomes 
>T ., (^ 



tf" 1 ) sin * - f (#' + #-*) sin 3* + ? (tf + #-) sin5*- .. . 

...... (167). 

Writing u=K'm. (163), ** = %v, and (165) becomes 



...... (168). 

(ii) Use of hyperbolic recurrence formulae. 

Write sin ff = sn (ft A') = tanh <&,', so that (159) gives 

n = - V sn a (iK' - *) = cot a ff = l/sinh 11 ^', < ft < K ... (169). 
Commencing with & in (169), use the recurrence formula 

tanh (<*' -*') = (&/O tanh *' ............ (170). 

We then have 

j8 = <k,'/2 n = *', sn 08, XT) = tanh ^ ', 
en 08, K) = sech <&,', dn 08, k") = sec 

and - i^ ftfc *) = 2 C B sin 

Equation (160) then gives 



As before, M is calculated from the recurrence formula 

tan (0 W+1 - <) = (^ n /) tan </>, 
giving a n u = <^ n /2 n = *. Equation (163) then gives for * ; the expression 

, _ sinh & cos^-g 2 sinhS^cosS^ + g 6 sinh 5$' cos 5<E -... 
cot v - c ^ n "^' 8 i n ^> __ ^2 cosh 3*' sin 3$ +V cosh 5* r sin 5O - ... 



Writing u = Kin (163), ** = Jir, and (171) becomes 
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SECTION X 

NOTE ON THE CALCULATION OF THE THIRD ELLIPTIC 
INTEGRAL IN TERMS OP THE COMPLEMENTARY A.G.M. 
SCALE 

When k is very nearly unity, it may in some cases be advantageous 
to make use of the more convergent complementary scale (<>'= 1* b '= K). 
The problem is to express the various functions appearing in the formula 



......... (174) 

in terms of recurrence formulae involving the complementary A.G.M. 
scale. The four cases discussed in Section IX, depending on the values 
of n, require separate treatment. The reader will have no difficulty in 
adapting the formulae of Section IX to the complementary scale should 
the occasion arise. We may note, in passing, the formula 

0(a- u) _ 



which will be found to be useful. 



SUMMARY OF FORMULAE 

NUMERICAL CALCULATION OF ELLIPTIC FUNCTIONS 
AND INTEGRALS 

[Sections containing new formulae are denoted by an asterisk *.] 
First Elliptic Integral. U = 

Second Elliptic Integral. E(j>, k) = J(l -t*sin> <f>)d<t>, E= 

Third Elliptic Integral. 
Jacobian Functions. 

sin </> = sn (u, k), cos < = en (w, k), J(l-te* sin 2 <) = dn (u, k). 

Scale of Arithmetico-Geometrical Means. 



1. Landens Scale of Increasing Amplitudes 

tan (<t> n+1 - </> n ) = (6 n /a n ) tan <f> n . 



(^"- E)IK= 
Z (+, k) = E(+, t) - 



G! sin <^i + c sin <^a + ... + c sin ^ B 
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*2. Calculation qf sn (u, k), cn (u, k), dn (u, k) in terms of argument u. 
Calculate $ n = 2 n a n u, n being chosen so that c n is negligible. Obtain 

^n-i"'&> &> ^o from recurrence formula sin(2</> n - 1 -^ w ) = - JL si 

a n 

Then sn (u, k) = sin < , cn (u, k) = cos ^ , dn (u, k) = cos ^ sec (<& - 

*3. 7%0 Hyperbolic Scale of Increasing Amplitudes 

tanh (<^ n+ i - n ) = (*n/n) tanh <^. Start with sinh ^ = tan <f>. 
'= Lt 



*4. Cakulation of sn (w, A'), cn (w, *'), dn (u, k') in terms oj argument u. 

(l/i) sn (iu, k) = sn (u, k')/cn (u, K), cn (iu, k) = 1/cn (u, k'), 
dn (iu 9 k) = dn (t^, #)/cn (w, A'). 

Calculate ^ n = 2 n a n u, n being chosen so that c n is negligible. Obtain 
<_! ... <^ 2 i ^u ^o from recurrence formula 

sinh (2^-! - < n ) = (c n /a n ) sinh <^ n . 
Then sn (u, k') = tanh <<>> cn (tt, *') = sech < , dn (u, k') = sech (& - < ). 



5. Calculation oj K'> E'> and the nome q = e- 
* * Hyperbolic recurrence formula tanh Xn+i 

= Lt [l/a*. X /2 n ], tanh Xo = N/F sin ^. 



. (Legendre.) 



log 2 = 2 (!)-> log -- log . (Gauss.) 

i a n+1 Ci 



a n 



6. Landen's Scale of Decreeing Amplitudes 



k') = log tan ft*- 
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x log tan ( JT 

<*n 

- ft) tan (2ft - ft) 



.... 

COS (2^+3-^+2) 

*7. Calculation of sn (#, '), en (w, '), dn (w, A') / terms of argument u. 
Calculate fti from sin ft t = tanh (a n w), T& being such that c n is negligible. 

Obtain ft^ ... ft, ^u *Ao from tan (^/i-i -^n) = -~ tan \ff n . Then 

7l 

sn (u, k') = sin ^ , cn (w, A;') = cos ^ , dn (t*, X:') = cos (2^ ^ ). 
*8. Hyperbolic Scale of Decreasing Amplitudes 

sinh (2i/r w+ i - W = ~ sinh ^. 
a n 

, k) - lja n tan" 1 (sinh i^ rt ), tan ^ = sinh ^ , 



i *) - ^(^ *) = tanh iAo cosh (2^ - ^ ) ~ sinh ^ - S' (An, c n ). 

jflL 2 

*9. Calculation of sn (w, X;), cn (u, k), dn (t^, X;). 

Calculate ty n from sinh ^ n = tan a n u, n being such that c n is negligible. 
Obtain ^ n _ 1 ... i/r 2 , ^, fa from taah (i/r n _ 1 - i/^ n ) =u- (c n /a n ) tanh i/^ rt . Then 
sn (w, ^) = tanh i^ , cn (^, >t) = sech ^ , dn (, k) = cosh (2^ - i^ )/cosh ^ . 

*10. Calculation of Third Elliptic Integral 4 cases. 

-rr / ? j\ sn (a ^') 

' 



^') f 

, . ; > v w^(a, X;) 

cn (a, k) dn (a, ) L 



Case I. - >fc 2 < w < 0. Case II. - GO < 2 < - L Hyperbolic Cases. 

'Method 1. (i) Circular, (ii) Hyperbolic 

recurrence formulae. 

Method 2. (i) Circular, (ii) Hyperbolic 
recurrence formulae. 



Case III. 1 < n < 
Case IV. < w < -f o> . 



tion. Case III, Method 1 (i). 
n = - #Ydn 2 (a, k') = - */(! - k' 2 sin 2 0), sin 6 = sn (a, #) < a < K 1 . 

1. Form A.G.M. scale (a = 1, 6 = '). 

2. Derive a w =- ^> w /2 n = * from tan (< n +i - ^ n ) = - n tan < n , commencing 

a n 

with <o = <^>. 

K 3 
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3. Derive a n a = logtan(^ + ^ n ) from sin(2^ n+1 -^ n ) = -r L si 

w 
commencing with fa = 0. 

oo 

4. Compute q from log q = 2 (^) n ~ 1 log (a n /a n+l ) - 

5. Compute ^ from 



_ 2# sin 2a w u sinh 2a n a + 2^* sin 4a rt w sinh 4a w a + 2q 9 + . . . 
~ 1 + 2^ cos 2a n tt cosh 2a w a + 2^ cos 4a w w cosh 4a n a + 2/ + . . . " 

(Legendre.) 
6. Then, 

K sin ^ + 5 (^, O]. 



APPENDIX 

In the preceding sections only those formulae suitable for numerical 
computation have been given. Many interesting and elegant properties 
of the arithmetico-geometrical scales in conjunction with recurrence 
formulae not specially adapted to numerical calculation are given below 
as problems. In some cases the main steps of the proof are briefly in- 
dicated, and unless otherwise stated the final results are new. 

Properties of the A.Q.M. scales (a - 1, b = '). 

1. Prove that a n = a n+l + <? n+1 , b n = a n+l -c n+l . Hence derive the series 



6 = ci n - GI + c 2 + c s + . . . . Gauss, G. W. Bd. in. p. 376. 
2. Prove that 



and hence *Ja n = N/a - \/c 2 - N/c 4 - ..., 

^n = \/6o + \/^ - Vc] - . . . . Gauss, G. JF. Bd. in. p. 376. 

3. Prove that c 2 n + 2c\ 4 . 1 = 2(a w -f-a w+1 )c WH . 1 , and hence derive the 
expansion 

(K- E}IK= (a + i) c x + 4 (a 2 + a,) c 3 

-f 16(a 4 + a 5 )^ 5 -f ...-f2 2n (a 2w -fa 2n+1 )c 2w+1 + .... 

4. Prove that 2b n -iC n = ^c 2 n _! 2c 2 n , and hence derive the series 



5. From the results of Ex. 4, prove that 

6 n > 1 c n -f ...), 



and hence the identity 

l=<* n a + 2 
6. From the recurrence formula of Ex. 4 prove that 



and hence a n z = 6 2 + ^c 2 - (cj 2 -4- Ca 2 + c 8 a +.-)> 

an 2 = a 2 ~|Co 2 -!-fe 2 + c 2 2 + c 3 2 -f-...). 

Gauss, 6?. fT. Bd. in. p. 376. 
7. Derive the recurrence formula 

log (0n-iAO = log (6/6-i) - 2 log (6+i/6 n ) 
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and hence prove that 



- log (o/i) - i log (,/*) - ... 
-(l/2 w )log 
JT/ AT - log (46,/Co) + f log (M>0 + I log (W + . . 



Gauss, G. W. Bd. m. p. 377. 
[Note that \irK'IK= (l/2 n ) log (a n /Cn) as n increases indefinitely.] 

8. Derive the recurrence formula 

(l/2 w ->) log KK-i) = (3/2 n+1 ) log (a n /b n ) 

+ (l/2 +1 ) log (a^J - (1/2*) log (0,-A-i)- 
Hence prove that 

TT/T/Jlf = log (160<AAo 2 ) + 32 (1/2'Olog (ajbj. 

Jacobi, JP 7 . # 52, p. 20L 

9. Show that the A.G.M. scale (a = 1, b Q ~k') may be calculated from 
the trigonometrical recurrence formula 



commencing with sin - k, sin ^ ~ (1 - A')/(l + k'\ etc. 
Hence derive the following formulae : 
(i) c^Jit 



from which, in general, expressions for c n and Cn may be formed. 
[Note that c n /a n = sin fl n , 6 rt /a H ~ cos O n .~\ 



(ii) log (atf/ir) - log (1 + sin ) - 2 log sec \ 0^. 

i i 

(iii) log {ir/(2JT^)} - J 2 log sec ^.^ 

(iv) j irJT7JSr= - J log q = J log (4/sin ft) - 1 (i) w log (1 + sin n+l }. 
[To derive (iv), make use of formula (61).] 

10. Show that the A.G.M. scale (a = 1, &o = #) may be calculated from 
the hyperbolic scale of decreasing amplitudes 

tanh 2x+ 1 - tanh 2 x, or x*+i = i log cosh 2x n , 
commencing with tanh 2x = *, or Xo = i log {(1 + *)/(! - )}. 

Hence derive the following formulae : 

(i) o^^tanhxo, c 2 = \k tanh x tanh Xi, c 8 =|Atanhxotanhxitanhx 2 > 
c? = } & tanh 2xi , cf - T V^ tanh 2xi tanh 2\ 2 , 
c 2 = A- if tanh 2xi tanh 2x 2 tanh 2x 3 , 
with similar expressions for c n and c w 2 . 
[Note that 6 n /a n = sech2xn, 
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(ii) log (2K/7r) = log (1 + tanh 2 X n). 



(iv) ^K'/K= log g~ i = log (2 coth Xo ) - 2 Q) n log(l + tanh x +i). 

[To derive (iv) make use of formula (61).] 

11. Show that the A.G.M. scale (a = 1, 6 = *') ma y ^e calculated from 
the hyperbolic scale of increasing amplitudes 

tanh 2^ n = tanh 2 \l/ n +i 9 or i^ w = |- log cosh 2^ n +i> 
commencing with tanh 2^ = &', or i/^ = ] log {(1 + #)/(! -*')} 

Hence derive the following formulae : 

(i) c n+l /c n = JU- 2 **, 6/a H = tanh 2i/r n , c rt /, t - sech 2^ w . 

(ii) log (2K/v) = 2 log (I + sech 2/r n ). 



(iii) log>/(2lT* f )} - 2 log coth <A,m. 



(iv) JirlT7Jr=logg"^ = ^/2 n " 1 , as ra tends to infinite values. 
[To derive (iv) make use of formula (60).] 

1 2. Writing = J { 1 + ( 1 - e ~ 8 ^n)4^ show that the recurrence formula 
of Example 1 1 may be written 

2</r w+1 - 4$ H = log 2 + log ^. 
Hence, denoting j8 = (1 - \/F)/(l 4- \/A;'), derive the expansions 

(i) log g- 1 = log (2/jS) + i (^-) n+1 log 4 +1 

(ii) log? -lo g (^) + 2(^) 4 -fl3(^) s + ^a/3) 12 +.... 
(iii) g =^ + 2(j^ + 15(jj8y + 150(}+.... 

(See also Example 24.) 

13. Derive the following series for TT from the formulae of Example 12, 
7r = log^- 1 -2.r 4 - nsf-^aP- ..., where *= J (2*-l)/C2* + 1). 
[Write k = k' = 1/^/2, then JT- ^ and log g- 1 - TT. Note that 

ir- I = 2314086, log a?- 1 = 3141602..., 2^*='0000070, 
so that two terms of the above series give the value of TT correctly to 
the tenth decimal. The remarkable approximation represented by the 
first term of the series was first pointed out by Legendre, F. E. t. I.] 

14. Given the nome g, compute the moduli k, k'. 
[The recurrence formula of Example 11 may be written 

^ n - 2^ = - log 2 4- log (1 - 
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Note that when n is sufficiently large ($T* l $ M = ?\agr l - We then 
have 

4^ w = 2 n log q- 1 - log 2 - 4*"\ 

Choose n so that the last term is negligible compared to the first two. 
The repeated application of the above recurrence formula enables us 
to calculate successively ^ n -i ... <A 2 and finally ^ in terms of which 
K = tanh 2^ , k = sech 2i/v] 

15. Derive the differential relations 

1 1 



for all values of n. Gauss, G. W. Bd. m. p. 380. 

16. From the results of Example 15 obtain directly the formula 



[In terms of the A.G.M. scale ( , b 0t c ) and the complementary scale 
(</, *o' ^o')> the formulae of Example 15 give 



But - ', ^ = c ', c = #0' and 2 = i 2 + c 2 , from which it is easily 
proved that 



Hence we conclude that ^ d (^ log |- J + a c? (^ log p J = 0. 

We also have the identity , d ( * ) + - , rf ( - 11 ) = 0, from which it 

a n a \a n / a n * \a n j 

follows that log ~ = C n -, , O n being a constant independent of 

* V n Q>n 

(a ny b n , c n ), (a n ', fe n ', c n '). In terms of Gauss 1 notation for the A.G.M., 
we note from Section III that 

a n ' = M(a ot c ) = 2 W 1/ (a n , c n ) = 2 n a n M (1, c n /a n ), 
while a n = M(a Q , i ) = M(a n , bn). Thus when TI becomes indefinitely 

great, C n = j' ~ log g = ^(^ log | . But if we write *' = <,/* 

Jlf(l, c tt /a n ) = ^7r/JE' = ^7r/log(4/A') as # becomes small. Thus, in the 
limit as n increases indefinitely, n = -jr, so that 
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17. Making use of the result of Example 16, prove that each of the 
ratios in Example 15 is given by 

V= Lt 



18. Derive Legendre's relation E/K+ E'lK' - 1 = Jir/CJEX") directly 
from the properties of the A.Q.M. scales (a = 1 , b - k f ) and ( ' = 1 , 6 ' - *) 
[Differentiating the identity log (a n bn) - log (tfn-i^n-i) = log( n /*n) and 
making use of the results of Example 15, we have 

d log (a n bn) - d log (a n - A-i) = d log (a w /6 ft ) = 2 H (0 n 2 /c 2 ) d log (aA)- 
Writing down similar equations for -l, w-2, ... 3, 2, 1 in place of 
9i and adding, we find 

*log(oA)-rflogfoA) = ( 2c > 3 + 4c * 2+ - + 2'V)cb- | rflog(aA) (i) 
and 

d log ('&') - d log (,'$;) = (2 Cl ' J + 4c a ' s + . . . + 2O Co'' 2 d log (o,'/V) 

...... (ii). 

Since a = ' = 1 & n d i ' = c o w h av ^ 

d log (oo'ftoO - d log (a 6 ) = d log (c /6 ) - a 2 o" 2 d log ( /6 ). 

[Ex. 1 
When n becomes large, a n -b n , so that 

dlog(a n bn)-dlog (.' V) - 2rflog (a w K') - - (4/ir) a n a n 'c 

[Ex. 17.] 

Note also the identity c ~ a rf log (o/*o) + c o'~ 2 ^ log (oW) = 0. [Ex. 16.] 
Subtracting (ii) from (i) and dividing out by 



we find 

(4/7r)a n a w '-a 2 

... + 2 n c n 2 + ...) - (2d' 2 + 4c 2 /2 



...... (Hi). 

But we have (K- E)fK= \ (c 2 + 2d 2 + . . . + 2 n c w a +...), o = i'/JPi o = 1 
with similar expressions for the complementary functions in terms of 
accented symbols. Thus (iii) becomes 



or, noting that c 2 + c /2 = 1, we derive Legendre's relation in the familiar 
form 

EIK+ E'lK' - 1 = i/( JOT). 

The equivalent of the series (iii) is given by Gauss (G. W. Bd. in. 
p. 380), but it does not seem to have been hitherto pointed out that 
Legendre's relation follows immediately from it.] 
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19. Verify the following identities : 



= 2 (1 4- 2^ + 2^4- 2# 18 + ...) 2 , 



Gauss, Nachlass, G. W. Bd. m. p. 447. 

20. Prom the results of Example 19, prove that if we form the 
A.G.M. scale commencing with 



c = A (2/4 + 2a?*- 
the wth term is given by 

a n = A<1 4- 2,r n + 2# n 4 + . . .) 2 , 6 = A (1 - 2# + 2# M 4 - . . .) 2 , 

C H - A (2^ + 2ar n * + 2^V + . . .) 2 , 
where log x n = 2 W log a? and A = JfcT(a , 6 ). 

Gauss, Nachlass, G. W. Bd. in. p. 448. 

21. Prove that ,r as defined in Example 19 is identical with the 
nome q when (a ~l, b Q -k'}. 
[When n becomes large, 

- 1 log q-frK-l K = Lt 1 log (fijcj 

n -^oo & 

- Lt ll 



so that ;r = <7.] 

22. From the results of Examples 20 and 21, derive the following 
expressions : 

Jv = Lr_?_l + 2<r / 



[From Example 20 we obviously have a = 1, 6 = ', c = i, h = 

23. If we denote = (1 - V^)/(l + V#), derive the following ex- 
pressions : 

W - (q + * + ^ +-)/(l + 2^ 4- 2,/ 8 +...), 
^T= 2ir (1 4- 2tf* + 20 W + 
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24. From the preceding examples derive the well-known expansion 

... etc. (See also Example 12.) 



25. If we write 



where log# n = 2 n log#, prove that the recurrence formula 

tan ^ n+1 - (OH+J/A^) tan 2ft, 
leads to the result 



- 

where tan ft, = *Jk' tan <, and K is given by the ^-series of Example 22. 

26. In terms of Landen's scale of increasing amplitudes (Section IV) 
tan (4Wi - <) = (6/a n ) tan < w , derive the following expansions : 
dn (u, K) = A (< , ) = n + Ci cos <fo + c 2 cos < 2 + . - . + c n cos <^ n + . . , 

oo 

^ sn ucnu = 2 2 n ~ 1 c 2 n _ 1 sin < n cos 



27. When (c n+l /a n+l ) is sufficiently small, prove that <^ n+1 of Example 
26 may be computed from the expansion 



2 sin 4<#> w - i (c n+1 /a w+1 ) 3 sin 

28. Prove that the recurrence formula of Example 26 is equivalent 
to A n = J (AVj + b 2 n)/^ n - lt and hence show that 

log [ ()/ (0)] = i -i, log (a n /A n ). 

[Jacobi, Fundamenta Nova (1881), 52, p. 204.] 

29. Prove that the recurrence formula a n tan <f> n+1 = A n+1 tan < M , where 
A w = *J(dn cos a ^> n + 6 2 sin 2 <t> n \ may also be written in the forms 




30. Prove that as a result of the recurrence formula of Example 29, 
rf^ n+1 /A n+1 . Hence derive the following formulae : 

= Lt (< n /a), where <#> = <<>, sin ^ = sn (w, *), 

n-^oo 

= 2ci cos <#> sin <k + 4c 2 cos <#>i sin < 2 + + 2 n c 



log [ (n, *)/e (0, *)] = 22- log [2a w /( n + A n )]. 
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[The recurrence formula which lies at the basis of the results obtained 
in Examples 29 and 30 was originally given by Gauss in 1818 in his 
paper " Determinatio Attractions, etc." G. W. Bd. m. (1866), p. 353, 
and was further discussed by Jacobi, " Zur Numerischen Berechnung 
der Elliptischen Functionen," G. W. Bd. I. (1881), p. 351. It is easily 
seen that these formulae are not adapted to rapid numerical calculation.] 

31. If sin <t> = sn (u, fc) prove that 

AI&Z ... A n ... , , 
tan a n u = - tan <<> 



[The result follows from Gauss' recurrence formula of Example 29 ; 
see also Jacobi, F. N., 36, G. W. Bd. i. p. 154.] 

32. Prove that Landen's recurrence formula 

tan fy n - </r n+1 ) = (<WiK+i) tan fa +1 
may be written tan 2\n+i = \/(^n/) tan \ n if we write 

tan \l/ n = 



33. If we make use of Landen's scale of decreasing amplitudes 

sin (Wn+i - fa) = (MO sin fa, 

commencing with fa = TT, prove that K'= (2K/Tr) log tan (TT + J^ n ) when 
n increases indefinitely. Hence show that the nome q may be computed 
from the formula q = cot 2 (JTT -f $fa) = (1 sin <A n )/(l + s ^ n 'An). Show, 
furthermore, that Legendre's relation leads to the identity 
a n (1 - sin ./O - (\*IK) {1 - tanh 

where 2 (^ n , c n ) is the series of equation (75). 



34. Prove that if the A.Q.M. scale (# , # ) is computed, the circum- 
ference / of an ellipse of semi-axes (#, ft ) is given by 

/ = (27r/a n ) (a, 2 - 2c 2 2 - 4c 3 2 - 8c 4 2 -...). 

35. If we denote by TI and r 2 the greatest and least distances re- 
spectively of a point on one circle to the circumference of a coaxial one, 
show that in terms of the A.G.M. scale (a = r lt 6 = ^a) the coefficient of 
mutual induction is given by 

M = {2**/a n } [c? -f 2c 2 2 + 4c 8 2 4- . . . + 2 n ~V +...], 

or, when the circles are very close together, in terms of the comple- 
mentary scale (ad = n + r a , b ' = r^ - r a ), by the formula 



[King, L. V., Proc. Hoy. Soc. A, Vol. c. 1921, p. 63 et seq. Several 
numerical illustrations are given in this paper.] 
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